arXiv:nucl-th/0503053v3 8 Aug 2005 


Electron scattering on halo nuclei 

C.A. Bertulanj] 

Department of Physics, University of Arizona, 1118 E fth St., Tucson AZ 85721 

Abstract 

The inelastic scattering of electrons on weakly-bound nuclei is studied with a simple model based 
on the long range behavior of the bound state wavefunction and on the effective-range expansion 
for the continuum wavefunctions. Three mechanisms have been considered: (a) dissociation of 
halo nuclei by high energy electrons, (b) dissociation by electrons present in a fixed target, and (c) 
Coulomb dissociation. It is shown that the properties of halo nuclei can be studied in electron- 
radioactive beam colliders using the electro-disintegration process. A comparison with fixed-target 
experiments is also performed. 
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I. INTRODUCTION 


A study of properties of weakly-bound neutron-rich, or halo nuclei, has been carried 


out intensively worldwide during the last decades |l|. Because of their short beta-decay 
lifetimes, halo nuclei are often studied in fragmentation facilities, where they are produced in¬ 
flight. The probes are hadronic, usually stable nuclear targets. Typically, one uses Coulomb 


dissociation, stripping, elastic scattering, etc. as nuclear structure probes. Such studies 
are complicated because the reaction mechanisms are not as well understood as with stable 
nuclear projectiles. The use of electromagnetic probes, e.g. electron scattering, is thus highly 
desirable. In fact, new experimental facilities for electron-scattering on unstable nuclear 
beams are under construction P| . An accurate determination of charge distributions in exotic 
nuclei can be obtained with electrons using inverse kinematics in a electron-nucleus collider 


mode 


Q. 


Electronic excitation, or dissociation, of nuclear beams can also be exploited for 
a deeper understanding of their structure. 

It is the aim of this work to explore basic results of electron scattering on the simplest of 
all nuclear halo structures, namely, a one-neutron halo system. The physics mechanisms and 
the conditions for the realization of electron scattering experiments are assessed. Such study 
has also an impact in nuclear astrophysics as it allows to deduce what are the lowest binding 
energies of halo nuclei possible in stellar environments, where free electrons are available. 

A high energy beam of weakly-bound neutron-rich nuclei dissociates as it penetrates a 
target due to the interaction with the atomic electrons. Since a heavy element target, e.g., 
^°®Pb, contains almost 100 electrons, the dissociation cross sections are large, assuming that 
each electron in the atom scatters independently on the projectile. Moreover, due to the 
atomic orbital motion, the innermost electrons have large relative energy with the incom¬ 
ing nucleus, increasing the dissociation probability. This process is of crucial importance 
in designing experiments aiming at studying properties of halo nuclei with the Coulomb 
dissociation method. 

The dissociation of neutron-rich nuclei, with small neutron separation energies, in stars 
can impose stringent limits on the stellar scenario where these nuclei play a role. For 
example, if the r-process proceeds partially out of equilibrium, the neutron radiative capture 
cross sections would have to be large enough to match the electron dissociation cross sections, 
with the appropriate neutron and electron density weights. 
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II. ELECTRON SCATTERING ON NEUTRON HALO NUCLEI 


I will consider the process e+a —>■ e' + b+c at small momentum transfers, q = (p^—p) /h, 
such that qR 1 (i? is the nuclear size). For simplicity, particle b is taken as a neutron and 
c as a core (inert) nucleus. The results obtained here are general and can be easily extended 
to the case of two-neutron halos. 

The differential cross section for this process is given by jj] 


dcTp 


2e2 


p 


2Jf + 1 EE' -I- c^p ■ p' -f- rrigC^ 


\pfi (q)l^^/, 


( 1 ) 


{hcY \p J 2 Ji + I q 

where E{E') and p(pO are the initial (hnal) energy and momentum of the electron, respec¬ 
tively. Ji{Jf) is the initial (hnal) nuclear spin, and Nf is the density of hnal states of the 
nucleus. The nuclear form factor pfi (g) is given by 


Pfi (q) = / Pf^ (r) d^r , 


( 2 ) 


where pfi (r) = ipi is the nuclear charge transition density, with rpi (ppf) equal to the initial 
(hnal) nuclear wavefunction. The cross section given by eq. ^ only includes longitudinal 
(also called Coulomb) excitations, dominant at low energy transfers ^,0,0]. 

Eq. HI is based on the hrst Born approximation. It gives good results for light nuclei (e.g. 
^^C) and high-energy electrons. For large-Z nuclei the agreement with experiments is only 
of a qualitative nature. The ehects of the distortion of the electron waves have been studied 
by many authors (see, e.g. ref. 0,0). For a rough estimate of this ehect, I follow ref. jli] |. 
For transition densities peaked at the nuclear surface with radius Rq, the correction due to 
Coulomb distortion is approximately given by 

dcTBorn/dhl 1 
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with 


P = 


120 




160a:3 


3 10 

1 -I- - cos (2a;) -|- 3a; sin (2a;) -I-(4-1-5 cos (2a;)) -I-a;^ sin (2a;) 

2 3 3 


X 


9 

-cos (2a;) 

4 


X 


oU 


i X 

^ 40 


( 4 ) 


where x = pRo/h, Si is the sine integral, Si(a;) = dtsint/t, and n, (a;) is the spherical 
Bessel function of the second kind. The above result is valid for monopole (I = 0) transitions. 
Corresponding expressions for higher order transitions are found in ref. [l0|. 
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TableHlshows the correction due to Coulomb distortion, eq. 01 for ^^Li and targets and 
several electron kinetic energies K^. One sees that below = 100 MeV it is important to 
account for Coulomb distortion of the electronic waves. A radius Rq = 3.5 fm was assumed 
for both nuclei. The Coulomb distortion correction decreases approximately linearly with 


Rq. 

An additional correction, due to nuclear recoil changes eq. [T] by a factor /rec — 
1 + {2Ex/M(?) sin(6'/2), where is the excitation energy, M is the nuclear mass, and 9 is 

the electron scattering angle. For the dissociation of weakly-bound nuclei, E^ -C Mc^ and 
this correction is much less relevant than the distortion of the electronic waves. I will neglect 
the Coulomb distortion and recoil effects from here on, bearing in mind that they should be 
taken into account in a more precise calculation. 


foe [MeV] 

QC'Li) 

Q^c) 

0.1 

0.879 

0.784 

1 

0.880 

0.786 

10 

0.887 

0.797 

10^ 

0.949 

0.903 

10^ 

0.994 

0.989 


TABLE I: The Coulomb correction factor, eq. ©, for electron scattering on ^^Li and and for 
several kinetic energies, (in MeV). 

In a simplihed model for the halo nucleus the radial parts of the initial and hnal wave- 
functions are represented by single-particle states of the form 

Ui{r) = Ai h^(i7]r) , Uf(r) = cos((5q) jif{kr) - sin((5q) ni^{kr) (5) 

where rj is related to the neutron separation energy Sn = /2^. hi^{irir) represents the 

large distance behavior of the bound state wavefunction, fi is the reduced mass of the 
neutron -f core system and hk their relative momentum in the hnal state, fo., ji^, and nq 
are the spherical Hankel, Bessel, and Neumann functions, respectively. A, is the ground 
state asymptotic normalization coefficient, which includes the normalization of the neutron 
single-particle wavefunction, and a spectroscopic factor which accounts for the many-body 
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aspects. This single-particle picture has 


□een used previously to study Coulomb excitation 


of halo nuclei with success jllL Il2l. llUL ll4L ll5L llfij | . 

The constant Ai (with spectroscopic factor equal to the unity) is used to normalize the 
bound state wavefunction, and corrects for the nuclear interaction range, Tq. In the case of an 
s-wave ground state (/j = 0), one has Ai = exp {’qr^) ^Jq/2TT (1 -|- qr^). For weakly-bound 
nuclei, 1/rj ^ tq and Ai ~ The ground state wavefunction entering the transition 

density integral, eq. |2l is well represented by the Hankel function hi.{ir]r). Note that the 
wavefunctions in eq. Elare not orthonormal. However, the transition density matrix elements 
of relevance for Coulomb excitation (and similarly for electron scattering) are dominated by 
the outside region (r > R) 13, Q|. Far from a resonance, the continuum wavefunction Uf{r) 
is small inside the nuclear radius. Its asymptotic dependence is well described by eq. El 

Using eq. El the form factor in eq. El can be calculated analytically by expanding 
into multipoles. The results will depend on the parameters S'„, R, and 5;^. To eliminate 
the dependence on R, the lower limit of the radial integral in eq. Elis extended to r = 0. 
The results for an s-wave ground state and the lowest order continuum angular momenta 
(// = 0,1, 2) are particularly simple. They are: 


(0) 
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where 


L = In 


yf + {k + qY 
R + {k- qY 


and M = tan 


-1 


k — q 


tan 


-1 


k + q 
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(7) 


In these equations = eZ(—1/H)^ is the neutron-core effective charge which depends 
on the transition multipolarity X {X = If for /j = 0). The effective range approximation 
^ 2 «+i has been used, where the parameters o/ and ri are the scattering 

length and the effective range, respectively. Notice that only for I = 0 the scattering length 
and effective range have dimensions of length. 


5 






















The I = 0 form factor has a large sensitivity to the orthogonality of the wavefunc- 
tions. If one assumes a zero-range potential for the neutron-core interaction, the scat¬ 
tering wavefunction, orthogonal to the bound-state wavefunction, is given by '0k^^(r) = 
exp(ik ■ r) — exp(z/cr)/[(77 -|- ik)r]. The s-wave scattering length is then just Oq = I/t], Using 
this value, together with rg = 0, in the equation for (g) leads to a large cancelation 
between the first and second terms. The / = 1, 2 form factors are also very sensitive to the 
scattering lengths and effective ranges. For example using oi ~ 5 fm^ and ri = 0 fm“^ 
reduces the magnitude of (q) by 10% for scattering at forward angles. These results 
show that it is very important to include the correct energy dependence of the phase-shifts 
to obtain an accurate description of electron scattering off halo nuclei . 

In what follows, I will use = e, i? = 0, /i = rriN (nucleon mass), and neglect the terms 
containing the effective-range expansion parameters in eq. El These approximations are not 
necessary but, with these choices, the numerical results will not depend on the charges and 
mass parameters of a particular nucleus; only on its neutron separation energy Sn- 

The total electron-disintegration cross section is obtained from eq. Q with the density 
of states given by Nf = d?k/ (27r)^, and integrating over k and Q. Figure Q shows the 
electro-dissociation cross sections obtained by a numerical integration of eq. ^ as a function 
of the separation energy, Sn, for electron bombarding energies equal to 0.1 and 0.5 MeV, 
respectively. One observes that only for very low neutron separation energies (5„ < 50 
keV) the electro-disintegration cross section becomes larger than 1 mb. If a more realistic 
model for I = 1 transitions is used, the cross section will be further reduced by: (a) a factor 
[ZajAaf' ^1/4 due to the effective charge, (b) by properly orthogonalized wavefunctions, 
and (c) by the energy dependence of the phase-shifts. 

Figure^also shows that the electron dissociation cross section increases appreciably with 
the electron energy. It is thus instructive to study the dependence of the cross section on the 
electron energy at high energies. The electron energy will be considered to be much larger 
than the energy transfer in the dissociation, i.e. E S> AU = E^ {E^ denotes the excitation 
energy). The scattering is peaked at forward angles and, from kinematics, q = k' cos 9 — k 
Ak ~ Ex/he. For energy transfers E^ of the order of a few MeV, one also has q p, p. 
Using eqs. El one obtains for the leading multipolarity (/ = 1) 



47ri e A qk 


( 8 ) 
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FIG. 1: Cross sections for electron induced breakup as a function of the separation energy, Sn, and 
for electron bombarding energies equal to 0.1 (dashed) and 0.5 MeV (solid). 


Using eqs. |H1 and [T] one obtains 
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48 ^ 2 ' 
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dQdEr 
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(9) 


The solid scattering angle can be related to the momentum transfer by means of dVl = 
2Tif^qdq/p^. The minimum momentum transfer for an excitation energy E^ is given by 
gmin = A/c = E^/hc, so that the integral over the scattering angle yields 


2 


2 





PH 


( 10 ) 


dE^ /ic2 E^ \E^^ 

Eq. fTTH shows that, for large p, the energy spectrum in electro-disintegration depends weakly 
on the electron energy through the logarithm function. This means that there is no great 
advantage (in terms of number of events) in increasing the electron energy when Ep ^ 
rUeC^- From eq. El one also sees that the energy spectrum increases sharply starting at 
Ex = Sn, peaks at Fa- = 8S'„/5, and decreases with Ex^^^ at large energies. This is the same 


characteristic spectrum as found in Coulomb dissociation of halo nuclei 


li 
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The integral over the excitation energy gives, to leading order, 


(Je {p) = GirV^e^ 



1 

fiC^Sn 


In 



( 11 ) 


For stable nuclei, with Sn — few MeV, the electron-disintegration cross section is small. 
The dependence of eq. on the inverse of the separation energy is most important for 
loosely bound nuclei. Using Sn = 100 keV, Eg = pc = 10 MeV, = e, and /rc^ = 10 ^ 
MeV, equation yields 25 mb for the dissociation cross section by high energy electrons. 
Note that the above equations are valid only if 3> rrieC^. They show that the electro¬ 
disintegration cross section increases very slowly with the electron energy. In contrast, as 
shown in figure ^ at low electron energies the cross sections increase much faster with E^. 

The arguments used here are only valid for dissociation (breakup) experiments. In the 
case of electron excitation of bound states, the matrix elements can become large for small 
excitation energies and cases where there is a large overlap of the wavefunctions. Conse¬ 
quently, the cross section can be much higher when these conditions are met. 


III. DISSOCIATION OF HALO NUCLEI BEAMS ON A FIXED TARGET 


A. Dissociation by atomic electrons in the target 


I use the Thomas-Fermi model to describe the electronic distribution in an atom. This 
approximation is well known, being described in many textbooks (see, e.g., ref. In 

this model, the electron density as a function of the distance from the atomic nucleus with 
charge Ze is given by 

. , , , 


p{r) = 




m 






, where x = br, and b = 2 ( — ] (12) 


Trip 


The function <F (a:) is the solution of the Thomas-Fermi equation 

$ 3/2 


(13) 


dx'^ ^ 1/2 

Numerical solutions of this equation date back to refs. excellent approximation 


was found by Tietz j2l|: 

$ (x) = 


(1 -|- ax) 


2 ’ 


where 


a = 0.53625. 


(14) 
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The probability density (normalized to Z) to find an electron with momentnm p is given 


by 

V{p) = |T>(p)|% where (p) = j ^p{r). (15) 

The electronic density p (r) has to be Lorentz transformed to the frame of reference of the 
projectile nnclens. Assnming a straight-line projectile motion with impact parameter b from 
the atomic center, the transformed density is 


p' (^) = IP + 7 ^ 2 ;^ 


(16) 


1 /Q 

where 7 = (1 — n^/c^) ' is the Lorentz factor, and v is the projectile velocity. 


The Fonrier transform in eq. ITHl becomes 
1 


D' (p) = 


i2n 


, 3/2 


dV e'P" = 


1 1 


where 


i27rf^ ^ 

P ={VoPzh) , r' = (b, 7 z) , 


Sr' , (17) 


(18) 


with pt [pz) being the transverse (longitndinal) momentnm. 

Since p{r) is spherically symmetric, eq. El can be rewritten as 

(P) = J dr r sin (Pr) y/p{r) . (19) 

For an atom at rest, very few electrons have orbital kinetic energies larger than 100 keV. 
In the case of only 3% of the electrons (2 electrons!) have kinetic energies larger than 
that. Bnt in the reference frame of a 100 MeV/nncleon projectile, 50% of the electrons have 
energies greater than 100 keV. 

Assnming that each electron scatters independently, the total dissociation cross section 
by the target atomic electrons is given by 

/ roo POO 

d^P p'(p) cTe (p) = 2ti dpt Pt / dpz V'{pz,Pt) (Te (p) ■ (20) 

^0 J—CO 

The separation of the above integral into longitndinal and transverse momenta is convenient 
becanse only the longitndinal momentnm component of the electrons is relevant for the 
dissociation of the projectile. 

Fignre El shows the dissociation cross section for a halo nnclens, with separation energy 
Sn = 100 keV, incident on a Pb target as a fnnction of the bombarding energy. Althongh 
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E|ab [GeV/nucleon] 


FIG. 2: Cross section for the electro-dissociation of a neutron-halo nucleus impinging on a lead 
target, as a function of the bombarding energy. A separation energy equal to 100 keV was used. 

the cross sections are small for incident energies eqnal to a few hnndred MeV/nncleon, 
they increase drastically as the bombarding energy becomes close to 1 GeV/nncleon. At 10 
GeV/nncleon the dissociation cross section is of the order of 1 barn. 

Gomparing the above resnlts with those obtained in section 2, we notice that there are 
different energy scales for electro-disintegration on hxed targets (by atomic electrons) and 
on a collider-beam mode. This is dne to the Lorentz transformation and to the large density 
(compared to an electron beam) of electrons in a heavy atom. Thns, with beams of halo 
nnclei with a few GeV/nncleon one conld, in principle, perform similar stndies as with 
electron-radioactive beam colliders. The disadvantage is that the Goulomb dissociation 
cross sections of loosely-bonnd nnclei are mnch larger, as shown in the next section. 


B. Coulomb dissociation 


Gonlomb dissociation of halo nnclei has been considered long time ago 




For the 


leading electric dipole transitions from an s- to a p-wave, the Gonlomb dissociation cross 
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section is given by ^ [l^ 


da. 
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32 
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V^{E^-Sr, 
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. 3/2 


In 


■yhc 
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( 21 ) 


where Z is the nuclear target charge, 6 = 0.681..., and R is the strong interaction radius 
{R ~ Rp + Rt)- Note the similarity with eq. ^1 in the dependence on the excitation 
energy E^, because the dipole operator is the same in both cases. However, the argument 
of the logarithm is different because of the small electron mass. Moreover, the coherent 
electric held of the projectile yields a factor which substantially increases the Coulomb 
dissociation cross section for large-Z targets. 

The total cross section for Coulomb dissociation as a function of the bombarding energy 
(i.e., as a function of 7 ) is given by 


27r 




--efr 




■jhc 


Using the same values listed after eq. for 10 GeV/nucleon projectiles impinging 
on Pb targets, yields cross sections of approximately 24 barns. This is much larger than 
that due to the dissociation by electrons in the target. But the contribution of the later 
process comprises 5% of the total disintegration cross section, and should be considered in 
experimental analysis. 


IV. CONCLUSIONS 

In this article the inelastic scattering of electrons off halo nuclei was studied, with em¬ 
phasis on the energy dependence of the dissociation cross sections. It is shown that the 
cross sections for electro-dissociation of weakly-bound nuclei reach ten milibarns for 10 MeV 
electrons and increase logarithmically at higher energies. This means that extracting infor¬ 
mation about the continuum structure of weakly-bound nuclei (e.g. scattering lengths and 
effective ranges, as in eq. El) can only be done if the intensity of the radioactive beam is 
very large, or if the collider allows for a large number of sequential interactions between the 
electrons and the nuclei at different crossing points. This conclusion can be drawn from £g- 
ured where a steep decrease of the dissociation cross section with S'„ is seen. Halo breakup 
experiments (common in hxed-target radioactive beam facilities) are difficult to carry out 
in electron-radioactive beam colliders, but not impossible if Sn is small. 
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A new facility is under construction at the GSI/Darmstadt, Germany. Experiments 
in a collider mode are planned so that electron beams will cross radioactive beams with 
center-of-mass energies of 1.5 GeV, i.e. 0.5 GeV electrons impinging on a 740 MeV/A 


counter propagating ion 


E 


24l |. For light, neutron-rich, nuclei luminosities of 10^^/cm^.s 


are expected. The approximate eq. [TT] yields cross sections of the order of 1 mb for S'„ ~ 1 
MeV, what means an estimated 100 events/second. 

I have also shown that electrons present in a hxed nuclear target access similar scattering 
conditions as in an electron-radioactive beam collider. However, Goulomb excitation cross 
sections are much larger in the case of a heavy nuclear target. In view of the scientihc impact 
of an electron-radioactive beam facility these results are useful for guidance in planning 
future experiments. The role of electron (and photon) scattering on exotic nuclei in stellar 
environments is also of interest for stellar modeling and work in this direction is in progress. 
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